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Among the few methods which have been proposed to create small-scale swimmers, those relying 
on self-phoretic mechanisms present an interesting design challenge in that chemical gradients are 
required to generate net propulsion. Building on recent work, we propose that asymmetries in ge¬ 
ometry are sufficient to induce chemical gradients and swimming. We illustrate this idea using two 
different calculations. We first calculate exactly the self-propulsion speed of a system composed of 
two spheres of unequal sizes but identically chemically homogeneous. We then consider arbitrary, 
small-amplitude, shape deformations of a chemically-homogeneous sphere, and calculate asymptot¬ 
ically the self-propulsion velocity induced by the shape asymmetries. Our results demonstrate how 
geometric asymmetries can be tuned to induce large locomotion speeds without the need of chemical 
patterning. 


I. INTRODUCTION 
A. Background 

Achieving self-propulsion at the micro-scale is essen¬ 
tial to many biological organisms and functions, includ¬ 
ing migration, feeding or escaping aggressions mm, and 
reproductive success of other larger species (e.g. mam¬ 
mals 0 ). From an engineering point of view, the design 
of self-propelled systems or micro-/nano-robots may of¬ 
fer important opportunities in particular for biomedical 
applications, to perform controlled therapeutic or diag¬ 
nostic tasks [10. 

Many efforts dedicated to the design of such microscale 
artificial “swimmers” have been inspired by the biolog¬ 
ical world, where viscous locomotion is achieved in the 
absence of inertial forces [6], for example using rotation 
or waving of rigid or flexible filaments. Experimentally, 
three broad categories of synthetic swimmers have been 
realized so far: (i) rigid [IHS] or flexible helices [TOHI^ or 
planar filaments m, inspired by bacterial flagella, forced 
by an external magnetic field in order to achieve propul¬ 
sion (see also Ref. m for a variation using nearby sur¬ 
faces); (ii) rigid bodies moving under the action of an 
external standing-wave acoustic field [Til [T5] : and (hi) 
so-called phoretic swimmers [TB] . The first two critically 
rely on the existence of an outside forcing in order to 
move, which may not only limit their applications but 
also disqualify them from achieving force- and torque-free 
propulsion. In contrast, phoretic (or fuel-based) locomo¬ 
tion, which is the focus of the present study, relies solely 
on the interaction of a rigid body with the solute content 
in the surrounding fluid. 

The ability to generate an effective slip velocity along 
a solid boundary outside a thin interaction layer in re- 
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sponse to a local tangential solute gradient is at the heart 
of classical phoretic physics. It originates from local pres¬ 
sure imbalance which are a consequence of short-range 
solute-particle interactions m- Classically, this phoretic 
mobility is responsible for the migration of inert particles 
in externally-applied chemical gradients. Specifically, a 
particle with uniform local surface mobility M placed in 
a far-held chemical gradient G of a neutral solute expe¬ 
riences a distribution of slip velocity on its surface lead¬ 
ing to a global phoretic velocity U = —M.G. For other 
phoretic mechanisms such as diffusiophoresis of charged 
solutes or thermophoresis, the slip and phoretic velocities 
depend instead on the gradient of the logarithm of the 
concentration/temperature, and this linear relationship 
is only observed for sufficiently small gradients [T7] . The 
basic idea of autophoresis is to combine such phoretic 
mobility with a chemical surface activity. Using chemi¬ 
cal reactions catalyzed at its surface, a chemically-active 
particle becomes able to generate the tangential gradients 
necessary to achieve its own propulsion [T5| . 

The feasibility of self-diffusiophoresis was recently 
demonstrated in several experimental studies using 
the catalytic decomposition of hydrogen peroxide on 
platinum-coated surfaces [TM2T] . although the exact 
physico-chemical mechanism at play is still under debate 
E5] . Note that this mechanism shares several im¬ 
portant similarities with self-thermophoresis |24j or the 
self-propulsion of droplets throngh Marangoni effects, for 
which the slip velocity is replaced with a surface shear 
stress discontinuity P5H?7] . 

In order to generate self-propulsion, breaking spatial 
symmetries is required. Indeed, the diffusion of a so¬ 
lute outside a homogeneous spherical particle leads to 
a spherically-symmetric concentration, and thus no slip 
velocity and no self-propulsion. Two routes have been 
identified to break symmetries and enable propulsion: (i) 
the chemical patterning on the surface of the particle, 
as used in most experiments with “Janus” particles [191- 
EHEH]; ( ii) an instability mechanism resulting from the 
nonlinear advective coupling of the solute to the phoretic 
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flows it creates, which spontaneously break symmetries 
and propel isotropic particles or droplets [271ES] ■ 

An alternative route to symmetry-breaking originates 
solely from geometry. Consider a particle with homo¬ 
geneous surface properties (i.e. uniform surface activity 
and mobility). In the absence of inertia, an asymmetry in 
the particle shape will in general create non-homogeneous 
solute (or reactant) concentrations along its boundary, 
which are then likely to have non-zero average and there¬ 
fore lead to propulsion. This idea is at the heart of exper¬ 
iments on collective phoretic dynamics |30j and it was re¬ 
cently the focus of an article analyzing the self-propulsion 
of a near-sphere with low-order azimuthal perturbations 
[31j using the osmotic framework of Brady and coworkers 
[32LI33j . In this paper, we tackle the same problem within 
the classical continuum framework of self-diffusiophoresis 
[34ljM] , focusing on two classical geometries amenable to 
analytical calculations. 


B. Intuitive model 


Before presenting detailed calculations in the following 
sections, we illustrate here intuitively the idea of acquir¬ 
ing locomotion from shape asymmetries by considering 
the case of a swimmer composed of two rigid spheres 
which share the same uniform chemical surface proper¬ 
ties. The spheres have radii i?i and i? 2 , and their centers 
are separated by a distance d. We are going to show that 
phoretic locomotion is guaranteed provided Ri 7 ^ i? 2 - 
Each sphere emits a solute with a uniform and identical 
rate A which diffuses in the fluid domain with diffusion 
constant k. We denoted by the unit vector joining 
the spheres’ centers (Fig. [^. In the limit of large sep¬ 
aration between the spheres, d ^ Ri, the leading order 
concentration field can be obtained by superposition of 
the distribution generated by each sphere independently 


c(r) 
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with corrections arising at higher order in 1/d. Each 
sphere is then exposed to two different phoretic contri¬ 
butions: firstly its own concentration field but since this 
is isotropic, it does not lead to any surface gradient or slip 
velocity; secondly the concentration gradients generated 
by the other sphere. Assuming a constant phoretic mo¬ 
bility Ai relating flow velocities to chemical gradients, 
the propulsion velocities of each spheres, U{ and U 2 , 
which would arise if they were individually force-free are 
given by 


f _ AMRl_ 
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A rigid two-sphere system (i.e. where d is kept constant) 
moves thus at speed U in a Newtonian fluid of viscosity 
rj such that the total hydrodynamic force is zero. In the 
far-field limit, the hydrodynamic resistance of sphere i is 


6TTfj,Ri and the total hydrodynamic force acting on the 
fluid is 67177 i? 2 (U — U 2 ) + 67177 i?i(U — U{) = 0 leading 
to locomotion at speed 

AAiRiR2{R2 — Rl) 

“ = ^d-{R, + R,) 

As long as 7 ^ R 2 , a net phoretic velocity is there¬ 
fore induced. The origin of this velocity is purely in 
the chemical asymmetries resulting from geometric dif¬ 
ferences, since both spheres are similarly chemically ho¬ 
mogeneous. The purpose of this paper is to illustrate and 
analyze these ideas further using two exact calculations. 


C. Outline of the paper 

In this paper, we study in detail how shape asymme¬ 
tries can lead to autophoretic locomotion. In order to 
combine physics with more practical considerations, we 
focus on two specific geometries. The first one, stud¬ 
ied in detail in (jllj determines the exact solution for the 
homogenous two-sphere system introduced above. This 
geometry allows us to address large shape asymmetries 
and can provide the basis for further experimental in¬ 
vestigations. The second geometry is that of a near 
sphere, addressed in m which allows us to quantify 
which asymmetric surface modes play a role in the transi¬ 
tion to locomotion and to address the question of optimal 
self-propulsion. 


D. Phoretic continuum framework and scaling 


Both problems in this paper are treated within the con¬ 
tinuum framework of autophoretic propulsion |34H36j , by 
considering an isolated rigid s yste m S (a two-sphere sys¬ 
tem in pi] or a near-sphere in (III I in an unbounded fluid 


domain of dynamic viscosity 77 and density p. The rigid 
body interacts with a solute species of local concentra¬ 
tion C that diffuses in the fluid medium with diffusivity 
K. The interaction layer thickness A is assumed to be 
small enough for the classical slip-velocity formulation to 
be valid m- The particle’s surface chemical properties 
are here characterized by a homogeneous surface activity 
A and mobility M.. At the surface of the particle, the 
solute is thus released {A > 0 ) or absorbed {A < 0 ) with 
a fixed flux, so that 


Kn-Vc = —A on 5, (4) 

and a slip velocity arises proportional to the local solute 
concentration gradient along the surface, following 

u = A4(l — nn) • Vc on 5. (5) 

If advection of the solute by the phoretic flows can be 
neglected (namely if the Peclet number, Pe = UTZ/k, is 
sufficiently small, with U the typical velocity magnitude 
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and R the characteristic size of the solid system), the 
solute has a purely diffusive behavior outside S, so that 
the concentration relative to the far-held solute content 
c = C — Coo satishes 

\7^c = 0. ( 6 ) 

Because the typical particle size and phoretic how mag¬ 
nitude are small, inertial effects in the how dynamics can 
be neglected, i.e. the Reynolds number Re = fUTZ/r] is 
small. The how held resulting from the phoretic slip at 
the boundary can then be solved for, in the reference 
frame attached to the particle, using Stokes’ equations 

? 7 V^u = Vp, V • u = 0, (7) 

together with the mobility condition in Eq. ([^. At in- 
hnity, we have 

u(r —>• oo) ~ —(U -I- r2 X x), (8) 

where (U, O) are the self-propulsion velocity and rotation 
of the solid system. As both the solid and huid inertia 
are negligible in the Stokes limit, the solid system must 
remain force-free, a condition that uniquely determines 
the self-propulsion kinematics (U, il). Since we focus in 
the following on axisymmetric systems, we have by sym¬ 
metry = 0 and U is parallel to the axis of symmetry. 

The problem is non-dimensionalized using R as charac¬ 
teristic length scale (the radius of one sphere in the two- 
sphere system or the mean radius of the near-sphere), 
\AM.\/n as characteristic velocity, \A\R/k as character¬ 
istic concentration fluctuation and r]\AAi\/TZK. as char¬ 
acteristic pressure. In the purely diffusive limit (Pe = 0), 
the only non-dimensional parameters arise from geome¬ 
try. 

This fixed-flux approach can be generalized to a simple 
one-step chemical reaction where the solute is consumed 
{A < 0 ) at the surface at a rate proportional to its con¬ 
centration m- In that case, A = —K,c where the reac¬ 
tion rate K, is now the uniform chemical property. The 
typical activity is now ICCao- This introduces an addi¬ 
tional non-dimensional parameter, the Damkohler num¬ 
ber Da = ICTZ/k,, which quantifies the limitation of the 
reaction rate by diffusion. For Da <C 1, diffusion is suffi¬ 
ciently fast to replenish the solute content in the vicinity 
of the sphere and the solute consumption occurs approx¬ 
imately at a fixed rate, while for Da ^ 1, the reaction 
is limited by the depletion of the solute content which 
cannot be compensated by the slow diffusion. 

II. AUTOPHORETIC LOCOMOTION OF A 
HOMOGENEOUS TWO-SPHERE SYSTEM 

A. Problem formulation 

We consider a system consisting of two spheres, iSi and 
S 2 , of respective radii i?i and i? 2 , with a fixed distance 


d maintained between their respective centers (either 
through long-range interactions or through a connect¬ 
ing rod with negligible hydrodynamic influence). In con¬ 
trast with existing studies focusing on a geometrically- 
symmetric dimer with chemical asymmetry |381139j , we 
specifically investigate here the effect of the geomet¬ 
ric asymmetry (i.e. different radii) for a chemically- 
symmetric system. Both spheres have uniform surface 
properties (activity and mobility) so that the axis joining 
their centers is an axis of symmetry for the problem. 
We first seek a solution of the diffusion problem for the 
solute concentration (relative to its far-held value) 


V^c = 0 outside the spheres. 

(9) 

n • Vc = —A on 5i and 52 , 

( 10 ) 

c(r —>■ 00 ) —?► 0 , 

( 11 ) 


with A = A/|A| = ±1, the dimensionless activity. Con¬ 
centration gradients at the surface of the two spheres 
generate surface slip velocities, and the following hydro¬ 
dynamics problem must then be solved 

V^u = Vp outside the spheres, (12) 

u = Ug = M(1 — nn) • Vc on and ^ 2 , (13) 

u(r —>• 00 ) ^ —Ue^, (14) 

where the swimming velocity U is such that the total 
hydrodynamic force on the two-sphere system is zero and 
M = \ = ±1 is the dimensionless mobility. 

Alternatively, using the reciprocal theorem for Stokes 
hows [40) . the swimming velocity can be determined di¬ 
rectly from the slip velocity distribution Ug on the two 
spheres as 

U = --^f Ug-cr*-ndS', (15) 

^ JSi,S2 

where (u*, cr*) is the solution of the dual hydrodynamic 
problem obtained by imposing a steady velocity 17* to 


the equivalent rigid (no-slip) two-sphere system 

V^u* = Vp* outside the spheres, (16) 

u = 17*02 on and ^ 2 , (17) 

u*(r —>• 00 ) —>■ 0 , (18) 

corresponding thus to a total hydrodynamic force F* 

F* = [ e^-cT*-ndS, (19) 

J (Si ,52 

with cr* = -p*l -I- (Vu* -I- Vu*"^) . 


B. Computing the swimming velocity 

1. Bi-spherical geometry 

Taking advantage of the geometry of the problem and 
its symmetries, a bi-spherical polar coordinate system 


4 



FIG. 1: Left: Notation for a two-sphere auto-phoretic system. The surface activity is uniform on both spheres, corresponding 
to a rate of solute release/absorption equal to A. Right: Bi-spherical coordinates: surfaces of constant r (solid) and /r (dashed) 
are shown; the surfaces r = r± correspond to the boundaries of the spheres. 


(r, /r, (j)) is used. Noting {p, </>, z) the classical cylindrical 
polar coordinates, the bi-spherical coordinates (—00 < 
T < 00 , —1 < ^ < 1) are defined as 


P = 


a^/l - 
cosh T — p' 


a sinh r 
cosh T — p 


( 20 ) 


In this system of coordinates, r = tq is a sphere centered 
on the z-axis at z = o coth tq with radius a /1 sinh tq |. The 
origin of the system of coordinates as well as the constant 
a > 0 are chosen such that r = r+ > 0 (resp. r = t_ < 0) 
corresponds to the surface of (resp. ^ 2 ) located in the 
upper (resp. lower) half-plane. The physical parameters 
(d, i?i, R 2 ) can be expressed in terms of (o, r_, t+) as 


i?2 = - 


sinh r_ 


= 


sinh 


d = a(cothT_|_ — cothT_), 

( 21 ) 


and provided that d > i?i -|- i? 2 , the solution of the in¬ 
verse system is unique. The unit vectors e,- and e^, re¬ 
spectively normal to surfaces of constant r and p, are 
defined by dx/dr = and dx/dp = with hr 

and the metric coefficients 


2. Solving the solute diffusion equation 


Within this system of coordinates, the Laplace equa¬ 
tion for the concentration field becomes 

d f 1 ac\ d f l-p^ dc\ 
dr \cosh T — pdr) dp \cosh t — pdp) 

whose general solution vanishing at infinity is im 

00 

c{t,p) = y/coshr - c„(r)£„(/x), (26) 

n =0 


with Ln the Legendre polynomial of degree n, and 


c„(t) 


a„exp 



{'T - T+) 


+bnexp 




(27) 


The constants, a„ and are determined by imposing 
the flux boundary condition on each sphere (noting that 
for T = r±, n(r±) = 


hr — T 1 hu — 

cosh r — 

Consequently we have 


(cosh T — p) — p^ 


1 — p cosh T v/l — sinh r 

^ ©2 — Sp ^ 

cosh T — p cosh T — p 

yjl — p^ sinh T 1 — M cosh r 

r “Gj. -|- • — 'Bp, 

cosh T — p cosh T — p 

and (e,-, 0 ^, 60 ) form a direct orthonormal basis. 


( 22 ) 

(23) 

(24) 


1 Oc 

Br ■ Vc at r = r±. (28) 

hr OT 


Substituting Eqs. (22) and ( [26| , and projecting the 
previous equation along the Legendre polynomial Lp{p), 
and using formulae summarized in Appendixj^we obtain 

coshr±c;(r±) - ^^c;_i(r±) - ^^c;+i(t±) 

+ sin^ ^^ Cp(T±) = ±aAv^e-(P+i/2)|r±|^ (39) 
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Together with Eq. (27), the previous equation form a 


linear system for the coefficients a„, which can be 
solved numerically after truncating the sum in Eq. (26). 


Once these constants have been determined, the slip 
velocity can be computed on the surface of each sphere 
as 

= —s/l- M2(coshr± - 


My'l- 




n—0 


(coshT± - 


3. Solving the dual rigid-body problem 


T±- 

m 


The solution of the general Stokes flow problem for an 
axisymmetric problem in bi-spherical coordinates can be 
expressed in terms of the streamfunction '0 (t,/ i). 




Ur = (coshr — Li) 

Ofl 


Un= - 


(coshr — d'ljj 




,2 dr ■ 
(31) 

The solution of the dual problem required to solve 


Eq. (15), i.e. Stokes flow vanishing at infinity with im¬ 
posed velocity ?7*e^ on the spheres, is given by m 

^ = (coshr - ^(1 - l^^)Lnil^)Un{T), 


n—1 


with 


(32) 


= 5^(i-m")l;(m)x 


(1 — /r^) sinhT± 

^ycoshT±^^]l 

U^{t±){ coshT± - n)- 


(35b) 


Projecting the two previous equations onto Lp{^) and 
Lp{tJ,) respectively, and using the relations given in Ap¬ 
pendix we finally obtain 

-p{p+ l)coshr±17p(T±) -b ' ^'^ Up-i{T±) 


+ 


2{2p-l) 

(p+l)(p + 2)(2p + 5) 


2{2p + 3) 


p- 1 


Up-i-liT±) = fp{p), 

(36) 

p -I- 2 


coshT±17;(r±) - 


--sinh T±C/p(T±) = fp{p), 

(37) 


with 


fp{p) = V2 

(38) 

re“(P+3/2)P±l Q-(p-^/'^)\r±\'^ 

fp{p) = V2sinh(T±) 


PCOSh'r± (p-l/ 2 )|r.(,| 
2p-l 


2p + 3 


2p-l 


(39) 


Using Eqs. (32), (36) and (37) applied at t± provides 


independent sets of four linear equations for the integra¬ 
tion constants a„, /3„, 7 „ and which can be obtained 
numerically. 


Unir) = a„ cosh -I-0 r-I-/3„ sinh -I-0 T 

-|- 7 „ cosh T -I- (5„ sinh t(33) 

In the previous equation, the four sets of constants a„, 
Pm 7n and Sn are determined by imposing the no-slip 
boundary condition on both spheres u* = 11*6^ in the 
dual problem or equivalently 


dip* , 17*(1 — pcoshT±) 

^ (coshT±-p)3 ’ 
dtp* U*(l-p^)sinhr± 

dr (coshr± —p)3 


(34a) 

(34b) 


These equations can be rewritten using Eq. (32) as 
(1 — pcoshr±) 


^Jco^LT±^^~p 


= E 


(35a) 


-n{n + l)L„(p)(coshT± - p) -b -(1 - p^)L'„(p) 


4- Swimming velocity 

From the previous two sections, the solute concentra¬ 
tion distribution and the dual problem streamfunction 
are completely determined. The swimming velocity can 
now be computed using the reciprocal theorem, Eq. (151. 


Because the phoretic slip Uj is purely along and the 
normal unit vector to the sphere’s surface is ±e,-, we only 
need to compute F* and = e^- • Vu* • -b • Vu* • e.^ ■ 
Using Eq. (34), on the boundaries of the spheres, we 
obtain 


Qq-*p(r±,p ) 


' = E L'^{p)Sn - coshr -b 


sinh"* r 


n—1 


2(coshT — p) 
(40) 


with 


Sn = -(coshT± - p)^/^U"(t±) 

+ — [sinhT±Un{T±) -b 3coshr±C/n(r±)] 

(41) 
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The total hydrodynamic force on the two-sphere sys¬ 
tem in the dual problem was computed in Ref. |41] as 

, A-Ky/2 ^ , 

F = -> n(n-|-l)(Q!„-|-7„), 

a ^ 

n—1 

so that the swimming velocity of the two-sphere 
is obtained as 

2^/ M^(T+,/r)g*^(T+,^)dA^ 

F* (coshr+-/r)2 

.1 (coshr_— J' 

In the previous equation, the integrals in /i are per¬ 
formed numerically (knowledge of the values of a„, 

Pm In and Sn Completely determines u® and a*^ on 
the boundary of the spheres). 



(42) 
system 

(43) 

(44) 


C. Results: self-propulsion of a two-sphere system 

The radius of the largest sphere is chosen as reference 
length scale, so that with no loss of generality we set 
i?i = 1 and i ?2 < 1- In the following, the self-propulsion 
properties and their dependence on d and i ?2 are exam¬ 
ined. When R 2 — 1, the system is up/down-symmetric 
and there is no net motion. Similarly, when R 2 1, the 
concentration distribution is only marginally impacted 
by the presence of the second sphere, and any net propul¬ 
sion velocity is infinitesimal. We thus expect the presence 
of an optimal ratio of sphere sizes. 


1. Far-field 

This is the result from the introduction, which we 
quote in dimensionless terms 

<F(R2 + 1} '■ 

When AM > 0 (particle with positive mobility releas¬ 
ing solute), in the far-held limit, the self-propulsion ve¬ 
locity is always oriented toward the smaller sphere. In 
this limit a maximum amplitude of the velocity is ob¬ 
tained for the optimal radius R 2 = — 1 « 0.41 and is 

equal to C/„iax = (72 - « O.lTAM/d^. 


2. Velocity for arbitrary distances 

When the size of the spheres is no longer small com¬ 
pared to d, higher order corrections in both the distri¬ 
bution of solute and the hydrodynamic held can become 
signihcant, and even dominant in the limit where the 
contact distance, dc = d — {Ri R 2 ) becomes small. 




FIG. 2: Top: Dependence of the self-propulsion velocity with 
the contact distance between the two spheres, dc, for R 2 /R 1 = 
0.75 (light red) and R 2 /R 1 ~ 0.5 (dark blue). Bottom: Same 
quantities plotted in log-log scale to show the scaling and 
compare with the far-held predictions (dashed lines). Here, 
both activity A = M = -\-l. 


In Fig. we plot the dependence of the swimming 
speed, [/, with dc when R 2 /R 1 = 0.5 and R 2 /R 1 = 0.75. 
Strikingly, the variation is non-monotonous and the ve¬ 
locity even changes sign. Specihcally, at a small contact 
distance, the self-propulsion velocity is positive (larger 
sphere in front) while at larger distances (and in the 
far-held) the self-propulsion velocity is negative (smaller 
sphere in front). As a consequence there is a hnite con¬ 
tact distance dc for which U = 0 despite the asymmetry 
in the geometry of the system. 

At large distance, the self-propulsion decreases as d~^, 
a direct result of the decay of the concentration gradi¬ 
ent created by each sphere on the other. This scaling 
is consistent with the far-held analysis, and the numeri¬ 
cal result obtained from the complete calculation shows 
a good quantitative agreement with the far-held asymp- 
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FIG. 3: Dependence of the swimming velocity with the size 
ratio, R 2 /R 1 , and the closest distance, dc = d — {Ri + R 2 ), 
between the two spheres. The dashed line indicates the con¬ 
figurations leading to no self-propulsion ([/ = 0). 


totic solution obtained in the previous section, Eq. (45) 
(Fig[^ bottom). 


When the spheres are nearly in contact, dc —0, the 
self-propulsion velocity converges to a finite and well- 
defined value, emphasizing that the concentration dis¬ 
tribution and flow field within the narrowing gap only 
has marginal influence on the dynamics of the entire sys¬ 
tem. In that limit, the two-sphere system behaves like 
a single solid with asymmetric (snowman) shape. We 
note that obviously our results are only applicable while 
the thin-layer approximation at the core of the classi¬ 
cal continuum framework remains valid (i.e. sufficiently 
large values of dc). However, because of the limited ex¬ 
tent of the contact region and its orientation relative to 
the direction of motion, we expect the results to be only 
marginally modified by a correction taking into account 
the finite size of the interaction layer. 


These results can easily be generalized to arbitrary size 
ratio, R 2 /R 1 , as plotted in Fig. Regardless of the rel¬ 
ative radius of the two spheres, self-propulsion occurs 
in different directions in the limit of near-contact and 
larger distances, and thus there always exist a finite con¬ 
tact distance dc « 0.1-0.2 for which no propulsion is 
observed. As expected, and optimal propulsion velocity 
is obtained for intermediate values of the size ratio (typ¬ 
ically R 2 /R 1 ~ 0.3-0.4), whose dependence with dc is 
weak. 


Within the {dc, R 2 /Ri)-p^3.iie, two optimal configura¬ 
tions can be identified. The first one, which is the global 
maximum, corresponds to two spheres in contact with 
R 2 IR 1 = 0.35 resulting in a velocity of 17 = 0.011. The 
second corresponds to a finite distance dc = 1.15 and a 
size ratio R 2 /R 1 = 0.31, resulting in 17 = —0.0098. 


3. Solute distribution 


The results above show that the distance between the 
two spheres critically impacts the propulsion velocity, 
particularly in determining its sign. To gain a better 
understanding of this effect, we show in Figthe solute 
concentration distribution for R 2 /R 1 = 0.35 (maximum 
velocity) and increasing distance. When A > 0 (solute 
release at the surface of the spheres), the solute concen¬ 
tration is always greater between the spheres than on the 
outside due to confinement: there the diffusive flux of so¬ 
lute can only take place on a reduced set of spatial direc¬ 
tions, leading to an increase in the solute concentration 
and in the gradient. This effect is even more pronounced 
when the two spheres are in contact, leading to singular, 
but integrable, solute concentration gradients near the 
contact point. 

Because of this solute distribution, the slip velocity 
is always oriented in opposite directions on the two 
spheres, regardless of their distance: when AM > 0 (resp. 
AM < 0), the slip velocity on each sphere is oriented 
toward (resp. away from) the other sphere. The con¬ 
tributions of each sphere to the propulsion are therefore 
always of opposite sign. In other words, if the spheres 
were independent they would move in opposite direc¬ 
tions, or, perhaps more quantitatively, the contribution 
of both particles to Eq. (44) have opposite signs. 

In order to quantify this more precisely, we define t/J 
with j = 1 or j = 2, the velocity of the two-sphere system 
obtained when only sphere j has non zero mobility (the 
other sphere’s mobility is set to zero). By linearity, the 
real self-propulsion velocity is obtained as U = U\ -\-Ul, 


and can be seen as a measure of the contribution to 
the total self propulsion. We plot in Fig. the depen¬ 
dence of the magnitude of U\ and with dc (note that 
we always have Ul < 0 < Ul). It shows, indeed, that 
rather than a fundamental change in the contribution 
of each sphere to the propulsion, the relative variations 
(and the slower decay of the contribution of the smaller 
sphere) is responsible for the change of sign in propulsion 
velocity. 

More specifically, for small dc, the contribution of the 
largest sphere dominates, a consequence from the large 
concentration gradients generated on that sphere near 
the contact point by the presence of the smaller one (see 
Fig. i left). This effect, mainly due to confinement, re¬ 
duces rapidly as dc increases, and the surface concen¬ 
tration distribution on the larger sphere loses its strong 
asymmetry. In contrast, the asymmetry of the concen¬ 
tration distribution on the smaller sphere is maintained 
at larger distance (Fig. right). 


4- Reactive vs. diffusive effects 

The chemical reaction at the surface of the spheres ad¬ 
dressed so far is a simple fixed-flux release/absorption 












FIG. 4: Solute concentration around a two-sphere system for a sized ratio R 2 /R 1 = 0.35 and three contact distances: = 0 

(left), dc = 0.24 (center), and dc = 1.15 (right). 



FIG. 5: Dependence of the magnitude of the contribution of 
sphere 1 (solid) and 2 (dashed) to the self-propulsion veloc¬ 
ity (defined as the velocity of the two-sphere system if the 
mobility of the other sphere is zero) with dc. 


of solute. A more general one-step reaction can be con¬ 
sidered by assuming that the solute is consumed with a 
uniform reaction rate, so that the dimensionless bound¬ 
ary condition on each sphere is replaced by m 

(46) 


Here the Damkdhler number Da = JCajK is a measure 
of the relative magnitude of reaction and diffusion. Note 
that Da = 0 corresponds to the previous situation (with 
A = —1). A generalization of this one-step fixed-rate ap¬ 
proach to chemical release (Da < 0) would not be phys¬ 
ically relevant as the release rate would be proportional 
to the local concentration, leading to a local exponential 
increase of the concentration and thus no steady solution 
to the diffusive problem. We exclusively focus on Da > 0 
in the following. 

In Fig. we plot the evolution with Da of the swim¬ 
ming velocity for the two optimal configurations iden¬ 
tified previously. In both cases, reactive effects are 
observed to significantly reduce the magnitude of the 
propulsion velocity and U ^ Da“^ for Da 1. As 
discussed in Ref. m, for finite Da, the solute consump¬ 
tion is limited by the reduction of its local concentration, 
which tends to reduce concentration gradient and slip 
velocities. 

For the optimal configuration where both spheres are 
in contact, a surprising change in the sign of the veloc¬ 
ity is observed at finite values of Da. While the two- 
sphere system self-propels in the direction of the smaller 
sphere at small Da (fixed-flux absorption, remembering 
that this is equivalent to the results of Fig. when cor¬ 
recting for the change of sign of the activity), finite-Da 
effects lead to a propulsion velocity in the direction of 
the larger sphere. 

This modification of the propulsion properties of the 
system for near-contact configuration is confirmed on 
Fig. where we plot similar results to those from Fig. 


n • Vc =1-1- Dac. 
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to the typical length scale of the particle. This prob¬ 
lem is similar to that presented by Ref. [3T], but using a 
different calculation framework. 


A. Geometry definition 

We consider here an axisymmetric particle of uniform 
surface properties, whose surface can be described in 
spherical polar coordinates as i?(/r) = 1 -I- with 

/r = cos 0 where 9 is the polar angle. We investigate the 
near-sphere limit, namely e 1, and consider the general 
case of a one-step kinetic reaction at the surface so that 
the activity and mobility of the particle are described by 

n • Vc =1-1- Da c, u = (1 — nn) • Vc, (47) 


FIG. 6: Dependence of the self-propulsion velocity of the two- 
sphere system with Da for the two optimal configurations 
for Da =0: R 2 /R 1 = 0.35 and dc = 0 (dark blue), and 
R 2 /R 1 = 0.31 and dc = 1.15 (light red). The results of the 
previous section are recovered when Da = 0 provided that 
the velocity sign is reversed {A — —1 here). 


for Da = 1 (left) and Da = 10 (right). When dc = 0(1), 
the effect of Da is a global reduction in the propulsion 
velocity (note the change of scales between Figs. and 
[^. However, when the spheres are close to each other, 
or even touching, the direction of propulsion is reversed 
in comparison to the case Da = 0. This effect occurs 
for smaller Da when the contrast between the sizes of 
the spheres is smaller. As a result, the existence of a 
non-propelling configuration for each size ratio is lost at 
finite and large Da, and most of the two-sphere swimmers 
self-propel in the same direction (namely with the larger 
sphere at the front), the direction being dictated by the 
slip velocity on the smaller sphere when A < 0: when 
reactive effects are significant, i.e. for Da = 0(1), they 
limit the role of confinement in setting the concentra¬ 
tion level and concentration gradients between the two 
spheres, thereby significantly reducing the contribution 
of the slip velocity on the larger sphere to the global mo¬ 
tion. 


III. AUTOPHORETIC LOCOMOTION OF A 
HOMOGENEOUS NEAR-SPHERE 


The previous section focused on a specific geometry, a 
system made of two spheres, which could be solved ex¬ 
actly even for large geometric asymmetries. We now turn 
to a different limit in which we address arbitrary pertur¬ 
bations in the shape of a spherical colloidal particle. In 
order to be able to compute the influence of each surface 
mode on the locomotion, we have to assume that the am¬ 
plitudes of the shape perturbations are small compared 


at r = 1 -F e^, and n is the normal unit vector pointing 
into the fluid domain 


er + eC\/l- 

■y/l -he2(l - /i2)^'2 


(48) 


We seek a solution of the autophoretic propulsion prob¬ 
lem as a regular series expansions in £ 1 of the solute 

concentration and velocity fields, and aim to compute the 
leading order contribution to the propulsion velocity 


c = + ec^ e^(? -I-... (49a) 

u = u°-heu^-|-£^u^-I- ... (49b) 

U = UO-hEU^-he^U^-h... (49c) 


Using Eq. (48), the expansion can also be applied to the 
unit vector normal to the surface of the particle 


n = no-I-£ni-I-£^n 2 -I-..., (50) 


e^, (51a) 

?'(m)\/ 1 - (51b) 

r(M)^(i-M^) ,, , r. -7 ... . 

--- imi (m) V 1 - M^ee.(51c) 

The deformation of the particle radius, ‘^(/i), is pro¬ 
jected onto orthogonal Legendre polynomials, so we write 
iid) = Note that by definition of the mean 

particle’s radius, ^0 must be zero. Also, at leading order 
only corresponds to a translation of the particle, its 
shape remaining symmetric. Hence, we do not expect 
any contribution from that mode to self-propulsion, at 
least at leading order. 


with 

no = 
ni = 

n2 = 


B. Solute concentration problem 

At all orders in £, the solute concentration c-^(r, fi) sat¬ 
isfies Laplace’s equation so that, enforcing the far-field 
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Ri! R\ 



FIG. 7: Same as Figure[^for Da = 1 (left) and Da = 10 (right). Since A = —I, the swimming direction is reversed. 


decay of c we obtain 


Lp[fi) 




(52) 


P—0 


where are constants obtained from the boundary con- 
dition at the surface of the particle, n • Vc = 1 + Dac, 
which can be rewritten at r = 1 using domain perturba¬ 
tions as 


-Dac° = 1 


dr 

or 


d'^e 

dr'^ 


(53a) 

.0 o 0 a 0 

-I- ^'(1 - -h Da^ —(53b) 


dfi 


dr 


dd^ 

dr 


— Dac^ = 'C^(l —M^)^t-1" 


2s^dc^ ^ ^''^{1 — fi"^) dc^ ^d'^c^ 


d^jL 


dr 


dr^ 


O„0 02 0 


2 dr^ U dr 


^2 q3^0 


g" d^C° 
2 dr'^ 


,(53c) 


Finally, Eq. (531 leads to the following solution for the 
second-order correction 


c^(dm) = - 


1 


^ ^ ^mnp^m^nRpif^) (56) 


1-h Da (n-I-1-I-Da )rP+i 

m,n,p—0 ^ 


where the third-order tensor Emnp is defined as 


E. 


2p+l 


mnp 


LpC^rnn^l^ 


/-I 


Qmn — 


(2n -f l)(n + 1) + n{n -h 4)Da -h nDa 


21 


n -h 1 + Da 


(n — 3 — Da) 


(1 - 


2(ti -t- 1 -f Da) 
where we note Lp = Lp{p). 

C. Stokes flow and swimming problems 


(57) 
RmLn 

(58) 


where all quantities are evaluated at r = 1. 

At leading order, the solution to the spherical problem 
is trivially obtained for c° as the isotropic solution 


c° = - 


1 


(1 -I- Da)r 


(54) 


After substitution of this result into Eq. (53), the first 


order correction to the isotropic solution concentration is 
obtained as c), = 2A^p/{p + 1) or equivalently 




^pLp(p,) 


1 -I- Da (p + 1 -I- Da )rP+^ 

p=0 ' ' 


( 55 ) 


The mobility of the particle, Eq. (47), imposes a tan¬ 


gential forcing on the fluid outside the particle result¬ 
ing in a global flow field in the Stokes regime. Using 
Eqs. (49), this boundary condition expressed at r = H-e^ 


can be converted, at each order, into a condition on the 
flow velocity on the spherical boundary at r = 1 

,_ gcP 

u° = -My/l - p'^—eg, (59a) 

dfi 

1 7,^ A- 2 f .d^c°\ 

dr \ dfi drdp J 

gfS ,_ 

+M{1 - —Br - M^' (59b) 
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= -e 




dr 


2 dr'^ 




ee 


dfi 

+M^'2(1 - fi^) + y/l^^^es^ 

+M^^Vl - “ \/l - 

where all quantities are again to be computed at r = 1. 

We now have to solve for the flow outside the unit 
sphere subject, at each order, to the boundary conditions 
at r = 1 and to the constraint that the total force on the 
unit sphere is exactly zero. This is strictly equivalent to 
solving for the flow outside the force-free non-spherical 
particle since the flow is also force-free (V • cr = 0) and 

I <T-ndS'— j cr-ndS' = j V-erdl/— j \/-crdV = 0, (60) 

Jr=l Jr=l+i Jvi JV2 

where Vi (resp. V 2 ) is the domain located inside (resp. 
outside) the unit sphere and outside (resp. inside) the 
actual particle. 

Replacing the flow problem on the unit sphere is partic¬ 
ularly convenient as it gives access to two major analyti¬ 
cal tools of low-Re swimming problems: (i) the reciprocal 
theorem to determine the swimming velocity U from the 
flow velocity on the unit sphere Ug, Eqs. (59), and (ii) the 


squirmer framework that provides an analytic solution 
for the flow velocity everywhere by projecting boundary 
conditions onto orthogonal squirming modes [42j . In par¬ 
ticular, the reciprocal theorem applied to a unit sphere 
provides the swimming velocity at each order from the 
surface velocity distribution as [ID] 


0{e) the surface velocity is purely tangential. Using the 
squirmer model framework [42l |43j , the flow field sat¬ 
isfying Stokes’ equations and the force-free condition can 
be written as 




2 


n(n -1-1) r 


(63) 


with '4’n(j') the radial streamfunction for the nth squirm¬ 
ing mode 

and a)j, the amplitude of that mode, obtained from the 
surface velocity as 


I = \ J ■eg)dfi 

Cnn{n - l)(n -k 1) 


1-1 

M 


1 + Da (2n -k l)(n -k 1 -k Da) 


(65) 


The squirming modes intensities characterize, among 
other things, the swimming properties and the far-held 
velocity created by the particle; ai = U is the particle’s 
swimming velocity, and the hrst mode also include a po¬ 
tential source dipole contribution. The second mode is 
related to the slowest-decaying singularity in the far-held 
of a force-free and torque-free particle, namely a stresslet 
The result in Eq. (65) conhrms that there is no net 


u^' = -{<) = Y 


swimming motion at hrst order (a\ — 0), but it however 
shows that the dominant stresslet is 0{e) and dictated 
by the second Legendre mode in A pure mode-2 

shape change would then lead to no net propulsion, by 
symmetry, but to a non-zero force dipole which will im¬ 
pact the bulk stress [H] and hydrodynamic interactions 
in a suspension. 

The sign of dictates the type of stresslet. When 


/ I rne sign 01 U 2 tiiciares rne lype 01 sriessiei. vvnen 

_ ^2 _ Yfxj dfi. (61) a 2 > 0 this corresponds to puller systems where the 
-1 ^ ^ thrust center is located in front of the draff center, simi- 


In jlllB we obtained that was isotropic. Conse¬ 


quently, it does not create any how at leading order, as 
expected, so that u° = 0. At 0(e), the boundary condi¬ 
tion on simplihes into 






,dc^ 

dfi 


ee 






-,dc^ 

dr 


ee 


(62) 


M 


1 -k Da 


a/1 - ; 


E 

.p—O 


i-p 


p -k 1 -k Da 




ee- 


Erom Eq. (^61| and using Eqs. (54), (55) and (59), we 


obtain = 0; hence, there is no net motion of the parti¬ 
cle at 0(e). However, a non-zero 0(e) how held is created 
through phoretic effects, and should be computed ev¬ 
erywhere in order to obtain from Eq. (59). At order 


thrust center is located in front of the drag center, simi¬ 
larly to the flagellated alga Chlamydomonas; in contrast 
0)2 < d correspond to pushers where the position of drag 
and force centers are reversed, as is the case for most 
hagellated bacteria. Eor M > 0 (resp. M < 0), pro¬ 
late particles (^2 > 0) will act as pushers (resp. pullers) 
while oblate particles (^2 < 0) will act as pullers (resp. 
pushers). As the solute consumption is enhanced near 
the poles of the prolate spheroid, when M > 0 the slip 
velocity is oriented away from the pole and toward the 
equator (puller), while for M < 0, the slip velocity is 
oriented from the equator to the poles (pusher). Inter¬ 
estingly, this dominant flow field decays rapidly with Da, 
a direct consequence of the slower diffusion: the reaction 
kinetics are slowed near the surface due to the depleted 
solute concentration resulting in smaller concentration 
gradients and slip velocities. 
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From Eqs. (631 and (65), the velocity gradient d\x"/dr 


at the surface of the sphere is now obtained as 
M 


—— r = 1 ) =• 

dr ^ 1 + Da 


C„n(n - 1); 


n—0 




2 n/W 

n + 1 


-L'M 


( 66 ) 


Substitution into Eq. (59) leads, for the velocity on the 
boundary, to 


U 2 = 


M 


1 + Da 


(Ymner + ^mn \/l^ — , by definition of n we have 


those with only even Legendre modes in their shape func¬ 
tion ^). One may also notice that as expected and 
do not contribute at this order as they merely change the 
radius or the center position of the sphere, not modify¬ 
ing its isotropy. More surprising is the role played by the 
second mode of deformation ^2 that dictates the dom¬ 
inant stresslet at 0{e): when Da = 0 (fixed-flux emis¬ 
sion/absorption), such a deformation does not contribute 
to the leading order propulsion, while it does as soon as 
Da ^ 0. 

When n S> 1, we notice that an = 0{n), and therefore, 
the infinite sum in Eq. (69) only converges if = o(l/ri) 
when n ^ 1 , a condition which is satisfied for regular 
shapes. If 7 denotes the angle between n and e^, then 


m,n—0 


with 


(67) 


n + 1 + Da 

( 68 a) 




n -I- 1 -I- Da 
(n — l)(2n — 1 -|- Da) 
n -I- 1 -I- Da 

71 + 1 + Da 


p—0 


p -I-1 -I- Da 


( 68 b) 


From the previous equation, the reciprocal theorem, 
Eq. ( [M] ), can be used to compute (see Appendix 
and at leading order we finally obtain O(e^) locomotion 
as 

00 

U = a„(Da)^„^„+i-f O(e^), (69) 

n—0 


tan 7 = i?'(/i)\/l — 


(72) 


It can be shown that a condition for divergence of the 
sum in Eq. (69) is the existence of a non-integrable sin¬ 


gularity in tan 7 . This, however, only includes marginal 
cases where the perturbative framework fails in the de¬ 
velopment of n as a regular perturbation series. 

A final interesting observation arises regarding the sign 
of the propulsion speed at different values of Da. From 


Eq. (69), we observe that for all n, the product a„(Da ^ 
1) • ^ I) is always negative. In other words, for 

any shape (i.e. any coefficients ^„), the particle will swim 
in opposite direction in the diffusion-dominated regime 
(Da 3> 1) and reaction-dominated regime (Da <C 1). 


D. Optimal autophoretic near-sphere 


We now use our asymptotic result, Eq. (69), in order 


to compute the optimal way to distribute the surface per¬ 
turbation modes maximizing the magnitude of the swim¬ 
ming velocity. 


with 


a„(Da) = 


F„ = 


[{n — 2)(n -|- 2) — (n -I- 7)Da — 3Da^]r„ 

(1 -I- Da)(2 -f Da)(n -I- 1 -I- Da)(n -1-2-1- Da) ’ 

(70) 

2Mn(n -|- l)(n — I) 


(2n-f l)(2n-f 3) 


(71) 


After having accounted for the difference in dimensional 
reference velocity, Eq. (69) is exactly equivalent to the 
result recently derived in Ref. m using a different, col¬ 
loidal calculation framework. 

This final result shows that it is indeed possible to cre¬ 
ate self-propulsion through shape asymmetries. In the 
case of infinitesimal perturbations of a spherically ho¬ 
mogenous autophoretic particle, the resulting swimming 
velocity is quadratic in the perturbation amplitude. The 
form of the result, Eq. (69), is consistent with the fact 


1. Optimization framework 

Formally, the leading order swimming velocity, 
Eq. (69), can be expressed in terms of the vector ^ = 


s{in)n as a bilinear form 

C/ = |-K(Da)-e (73) 

with the upper-diagonal linear operator K defined as 
A"mn(Da) = a7i(Da )dm,n-t-i ■ 

Since the swimming velocity is quadratic in the ampli¬ 
tude of the perturbation, optimization requires to impose 
some kind of fixed norm to guarantee that we remain 
within the perturbative framework in which this result 
was obtained. Hence, we define an objective function J 


that is zero for front-back symmetric particles (i.e. 




(75) 
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with A a Lagrangian multiplier and H a bilinear, symmet¬ 
ric positive definite operator corresponding to the partic¬ 
ular norm chosen for Maximizing or minimizing U is 
therefore equivalent to seeking solutions of 

(K + K^).^ = H.e with = (76) 


The previous equation is effectively an eigenvalue prob¬ 
lem for (4, A) with a constraint on the norm of the eigen¬ 
vectors In practice, for a given norm operator H, we 
seek the solutions of Eq. (75) for a finite number N of 
modes in This leads to a set of N numerical solu¬ 
tions, from which the shape with maximum velocity is 
extracted. 

The choice of the norm used to define H significantly 
impacts, as expected, the convergence of the optimal so¬ 
lution when TV —>■ oo. For example, a simple L 2 -norm of 
the perturbation amplitude ^ (denoted or of local 

slope perturbation (i.e. angle between n and e^, denoted 
which led, respectively, to the following tensors 


= /\l - (78) 


do not guarantee the convergence of the result as TV —> 
oo. In fact, in both cases, for finite values of TV, the 
optimal shape is obtained as a sharp cusp near the pole 
leading to a dominance of the higher-order modes in ^ 
over the more regular ones. This can be understood as 
follows. Bounds on the r.m.s. perturbation amplitude or 
r.m.s. perturbation angle do not rule out the presence of 
a spike of infinitesimal thickness and unbounded height 
on the surface, that would lead to a diverging of t/^ (the 
regular perturbation expansion of the propulsion velocity 
would be invalid). 

In order to guarantee well-posedness, and consistency 
with the perturbation approach, we choose a norm based 
on the r.m.s. curvature perturbation k = k + where k 
is the local curvature defined as 


Ws 


= Atn, 


(79) 


with t the tangential unit vector in the (e^, ee)-plane and 
s the curvilinear coordinate along the surface. At leading 
order in e, the curvature perturbation is obtained as 


k = + (80) 


and the corresponding H is obtained as 

Hmn = J (81a) 

hj = fiL'j -I- [1 — j{j + 1)]T„, Vj. 

Hmn is now a full matrix whose coefficients can be com¬ 
puted analytically (see Appendix [b|) . 



FIG. 8: Small-amplitude perturbations of a sphere leading to 
maximum swimming velocity (e = 0.3 was chosen for plot¬ 
ting purposes). There are two solutions, one convex (solid 
line) and one concave (dash-dotted line), leading to the same 
propulsion velocity U ~ O.Ole^TW. The unit sphere is shown 
for reference (dotted line). 


Note that because the propulsion velocity is a 
quadratic form of the shape perturbation, the transfor¬ 
mation ^ ^ leaves the propulsion velocity unchanged. 

Each optimal swimming velocity therefore corresponds 
to two different shapes obtained for e and —e. Also, the 
change —>■ (— 1 )”^„ simply performs a symmetry of 
the particle shape, changing the sign of the swimming 
velocity but not its magnitude. 


2. Optimal swimming shape for fixed flux (Da — 0) 


We plot Fig. the two optimal shapes leading to max¬ 
imum propulsion velocity at Da = 0 (fixed-flux absorp¬ 
tion). For positive mobility {M = 1), both shapes swim 
to the right [U > 0) with a small speed, U « O.Ole^. 
Both shapes are characterized by a sharp corner with a 
finite angle 7 at the pole, defined in Eq. (721. This angle 


is proportional to £, a consequence of effectively imposing 
the curvature change to be 0 (e). 

The optimal shapes appear to be cusped. It is a well 
known fact that, for example, in electrostatics, cusps can 
lead to divergence of electric fields, a result true in gen¬ 
eral to solutions of Laplace equation near tips or wedges 
in the geometry |45]. We believe that the same effect 
is observed in our simulations where the maximum ve¬ 
locity, resulting from surface gradients of the solution to 
Laplace’s equation, corresponds to a kinked geometry. 
As for the two-sphere system, we note that the frame¬ 
work used here is valid provide that the typical radius of 
curvature near the regularized cusp remains compared to 









14 



FIG. 9: Dependence of the swimming velocity on Da for the 
optimal shapes shown in Fig. [^(i.e. optimal shape at Da = 0) 


the thickness of the interaction layer. 

In Fig. we show the dependence of the swimming 
velocity associated with this optimal shape on the value 
of Da. In general, reactive effects are observed to reduce 
the swimming velocity, a property that was also observed 
for spherical Janus particles m and remain valid here: 
reactive effects tend to reduce concentration contrasts as 
the chemical reaction is slowed down in regions already 
strongly impacted by the depletion in solute resulting 
from the reaction. More quantitatively, at large Da, one 
observes that for a fixed particle shape U ~ Da“^. How¬ 
ever, this decrease is not monotonous: in fact, a reversal 
is observed in the swimming direction as already pre¬ 
dicted in the previous section. As a consequence, for a 
finite value of Da (Da « 0.94 in this particular case), 
the asymmetry-driven self-propulsion vanishes. 

3. Optimal swimming shape for arbitrary Da 

The previous result clearly suggests that the optimal 
swimming shape can not be independent of Da since any 
shape will perform poorly for some range of Da = 0(1). 
Our optimization analysis can hence be extended to fi¬ 
nite values of the Damkohler number and the results are 
presented on Fig. where we plot the overall shape 
(top) and the angles at both poles (bottom). For large 
Da, one recovers the general decay of self-propulsion ve¬ 
locities due to reactive effects and C/opt Da“^. The 
optimal shape is rather insensitive to fluctuations in Da 
when Da < 0.4 and Da > 0.7. Within a narrow range 
of Da, a sharp transition takes place. While at lower 
Da, the optimal shape presents a protrusion located at 
its front, for greater Da the protrusion is in its back. 
This change in the optimal configuration is related to the 
change in self-propulsion direction for a particle shape 
with a protrusion observed in Fig. For small Da, it is 
advantageous to have a protrusion in the front, while at 
larger Da, the protrusion should be located in the back 



Da 

FIG. 10: Top: Dependence of the optimal swimming velocity 
and optimal shape with Da. A front-back reversal in the 
optimal shapes occurs near Da ~ 0.5. The direction of the 
swimming velocity is indicated for each displayed shape by 
a green arrow. Bottom: Variation with Da of the summit 
angles on the optimal shape. 7 + (solid) refers to the right 
pole and 7 _ (dash-dotted) refers to the left pole. 


of the particle. 


IV. CONCLUSION 

The results presented in this paper for two ideal¬ 
ized setups demonstrate that geometric asymmetries 
are sufficient to enable self-propulsion of a chemically- 
homogeneous system. As in the work of Ref. [3T], this 
identifies an alternative route to self-propulsion with¬ 
out exploiting chemical patterning of the particle which 
might otherwise be practically difficult to achieve and 
control. We note however that the propulsion velocities 
are in general smaller than what would be achieved for an 
equivalent Janus system exploiting chemical asymmetry. 

The origin of this geometric self-propulsion mechanism 
stems from a combination of enhanced concentration gra¬ 
dients and change in the particle surface exposed to these 
gradients. The former corresponds to an enhancement 
of the concentration level and gradients, and hence slip 
velocities, near a protrusion from a surface of constant 
curvature, and is a generic feature of computing local so¬ 
lutions to laplacian on wedge or tip singularities |35] . The 
latter can be best illustrated in the near-sphere case: the 
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extension of the solid particle away from its mean radius, 
exposes its surface to a larger range of the concentration 
distribution created by the particle, which modifies sur¬ 
face velocities. The dependence of the final results for the 
swimming velocity on the various surface modes, on the 
deformation amplitude, and on the relative importance of 
reactive effects show however that it is a nontrivial non¬ 
linear interplay between enhanced gradients and shape 
changes which results in swimming. 

Most of our results were obtained with a minimal ki¬ 
netic model for the reaction at the surface of the parti¬ 
cle, namely a fixed-flux release or absorption, but we also 
showed that introducing a more complex one-step reac¬ 
tion kinetics does not fundamentally change the results, 
in particular the main idea of self-propulsion using shape 
asymmetries. Our work emphasizes however that the pre¬ 
cise kinetics and the relative importance of reactive and 
diffusive effects may significantly influence the sign of the 
propulsion velocity, and additional work will be needed 
on more complex, and realistic, surface chemical condi¬ 
tions (see for example Ref. [21]). Furthermore, while we 
focused on the idea of self-propulsion in this work, the 
same idea will be useful as a method to apply forces on 
stationary bodies, leading therefore to potentially new 
microfluidic applications in pumping and mixing |46j . 
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Appendix A: Useful properties of the Legendre 
polynomials 


We make use of the following in the main text of the 
paper 


2p+l 


y ^(1 - /r^)L(,Lpd/r = 5 , 


2p+l 

2p{p +1) J_ 

r-i 

|XZ/yjZ/pd^ — 


(Al) 


(A2) 


2p+lf'^ P+ljr P 


1-1 


2p + 3 


2p+l 


J ^{1-p‘^)L'^Lpdp = 


(A3) 

{p + 2){p + 1) 

2p + 3 
-ffhi,.,,-., (A4) 




2p+l 

2p{p + 1)7- 


y ^ /r(i - p^)L'^L'^dp= 


2p + l 
2 

2p + l 
2 



Lpdfi 

^ cosh T± — p 
pLpdp 
cosh T± — pL 


+ 2p - 

= y2e-(P+i/2)|r±| 


= v/2 


( Z±lp-(P+3/2)|T±| 

^^ 2 ^ + 3 


_p-(p-l/2)l'r±| 


2p- 1 


(A7) 


2p+l (1 - /i^)Lpd/r 

2p{p -I- 1) cosh T± — fj. 


= V2 


3-(p-1/2)|t±| 


2p- 1 
-(p+3/2)|r±| 

2p 3 


(A8) 


Appendix B: Computing the second-order swimming 
velocity 


Applying the reciprocal theorem, Eq. (61), to the sec¬ 
ond order surface velocity, Eq. ( 68 ), one obtains 
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where E^np was defined in Eq. (57) and 

and Ip are integrals of the Legendre polynomials 
that can be computed using classical properties of such 
polynomials and those listed in appendix |A| 
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Also, using the definition of Emni and the properties 
of Legendre polynomials, we obtain 
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with 

/„ =(n + 1) ^2nDa^ + 2n{n + 4)Da + 2(n + l)(2n + 1)^ 
+ n(n + l)(n + 2)(n — 3 — Da), (B8) 

dn =n^2nDa^ + 2n{n + 4)Da + 2{n + l)(2n + 1)^ 

+ n(n — l)(n + l)(n — 3 — Da) (B9) 


with a„(Da) defined in Eq. (69). 


Appendix C: Curvature norm tensor 


Substituting these results into Eq. (Bl), we obtain 
in the form 
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Finally, 
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n =0 


Choosing the r.m.s. curvature perturbation to define 
the H operator in Eq. (75), the symmetric tensor H^n 


is defined in Eq. (81) and is symmetric. In particular. 


using classical properties of Legendre polynomials (see 
Appendix [A| , it can be computed exactly , when m < n, 


as 
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